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We consider the consequences for the relic neutrino abundance if extra neutrino interactions are 
allowed, e.g., the coupling of neutrinos to a light (compared to m„) boson. For a wide range of 
couplings not excluded by other considerations, the relic neutrinos would annihilate to bosons at 
late times, and thus make a negligible contribution to the matter density today. This mechanism 
evades the neutrino mass limits arising from large scale structure. 
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Introduction. — The discovery of neutrino oscilla- 
tions means that neutrinos have mass, which requires 
physics beyond the Standard Model. The solar and at- 
mospheric oscillation experiments have measured neu- 
trino mass-squared differences b~m\ x ~ 7 x 10~ 5 eV 2 and 
bm\ 2 ~ 2 x 10" 3 eV 2 Q, which implies lower limits on 
two neutrino masses of \J Sm^i and \J Sm^ 2 - Since these 
oscillations have been shown to be dominated by active- 
flavor neutrino oscillations, the three neutrino masses 
are connected and become degenerate in mass if any are 
larger than \J 8m\ 2 3- Thus, at the present sensitivity of 
m„ < 2.2 eV (at 95% CL) the upper limit on neutrino 
mass from tritium beta decay applies to each of the three 
mass eigenstates. KATRIN, a proposed next-generation 
tritium beta decay experiment, will have sensitivity down 
to m v ~ 0.2 eV [ij. New neutrinoless double beta decay 
experiments will have even greater sensitivity, but only 
if neutrinos are Majorana particles 

Neutrino mass can also be measured with cosmolog- 
ical data. When neutrinos are relativistic, they free 
stream out of density perturbations, reducing the growth 
of structure. This results in a suppression of the matter 
power spectrum on all scales below that of the horizon 
at the time the neutrinos became non-relativistic, after 
which they act like cold dark matter. The extent to which 
this lack of clustering affects the distribution of matter 
today depends on the ratio of the energy density of the 
non-clustering component (neutrinos) to the total den- 
sity of matter. The former is 

p v = Y,m v n v = 93 ^ eV p cr , (1) 

where p cr = 3Hq/8wG is the critical density associ- 
ated with a flat universe; the total density in matter is 
parametrized as f2 m p cr . Here h specifies the Hubble con- 
stant, H = 100/ikmsec -1 Mpc -1 . The equality on the 
right in Eq. assumes the standard cosmological abun- 
dance. Recall that, in the standard scenario, neutrinos 
couple to the rest of the cosmic plasma until the weak 
interactions freeze out at T ~ 1 MeV. After neutrinos 
freeze out, their abundance scales simply as a -3 where a 
is the cosmic scale factor. Thus, in the standard cosmol- 



ogy, there are roughly as many relic neutrinos today as 
photons in the cosmic microwave background (CMB). 

Limits from structure formation on the sum of neu- 
trino masses now range from 0.5 to 2 eV, with the spread 
largely due to different assumptions about the relative 
bias between the mass and galaxy distributions Bias 
is one important issue, but this will be circumvented with 
future weak lensing surveys, which will measure the mass 
distribution directly. Indeed, it has been shown that 
these observations should realistically be able to reach 
the scale \J (JmL by which the discovery of neutrino mass 
is guaranteed These mass constraints depend on as- 
suming the standard relic neutrino abundance. Big-bang 
nucleosynthesis (BBN) constraints, combined with neu- 
trino mixing data, no longer allow the possibility of a 
significantly increased n v due to a large lepton asymme- 
try |8j. Are there other ways to alter the relic neutrino 
abundance, and specifically to lower it? 

If neutrinos have extra interactions so that they remain 
in equilibrium until late times, they would freeze out 
when they are non-relativistic, in which case their final 
abundance would be suppressed by a factor cx e~ m f/ T f , 
We show that new neutrino couplings in the allowed 
range can lead to a vanishing relic neutrino density today, 
hiding the effects of neutrino masses from cosmological 
observations. This possibility is falsifiable both directly 
and with other experiments. 

Interaction model. — We consider the cosmological 
consequences of coupling neutrinos to each other with 
bosons, through tree level scalar or pseudoscalar cou- 
plings of the form 

C = h,T',\',o + !h,V,-',u,o + h.c, (2) 

as in Majoron-like models, for example. The field 4> is 
assumed to be massless (or light compared to m v ). Vi- 
able models of this type have been discussed in Ref. |IJ. 
Here we assume that there is just one new boson, and 
that these new couplings are unconnected to the mech- 
anism of neutrino mass generation. Even tiny couplings 
can cause profound effects, as we will show. 

The solar neutrino 0, El and meson decay 0] lim- 
its on these couplings are very weak, \g\ < 10 -2 (here 
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and below we do not distinguish g or h type couplings, 
nor neutrinos and antineutrinos) . Neutrinoless double 
beta decay limits g ee < 10 -4 , but the other couplings 
may be much larger. Supernova constraints may exclude 
a narrow range of couplings around g ~ 10~ 5 , but the 
boundaries are model dependent Scalar couplings 

could mediate long-range forces with possible cosmolog- 
ical consequences, while pseudoscalar couplings mediate 
spin-dependent long-range forces, which have no net ef- 
fect on an unpolarized medium [l3j |. Since these con- 
straints can be evaded, and since in our case <f> couples 
only to neutrinos, we do not consider them further. 

The 4> boson can be brought into thermal equilibrium 
through its coupling to the neutrinos, and the v — <j> sys- 
tem may stay in thermal contact until late times, through 
the processes v<f> <-> V(f> and v <-> vc$>. Most important 
though is vv «-► (f>4>, a process which depletes the total 
number of neutrinos. In the standard case, the neutrinos 
decouple from each other and the matter at T ~ 1 MeV, 
but interactions with <p may keep neutrinos in equilibrium 
until they are non-relativistic, T ~ 1 eV, when the in- 
verse process becomes kinematically prohibited. In order 
to accomplish this, g must be sufficiently large; we show 
below that this requires g > 10~ 5 , well within the al- 
lowed range. While neutrino decay requires off-diagonal 
couplings, the effects considered here can occur with ei- 
ther diagonal or off-diagonal couplings. If the couplings 
are this large, all relic neutrinos efficiently annihilate into 
bosons, leaving no relic neutrinos today, thereby hiding 
the cosmological effects of neutrino mass. 

Past models of invisible neutrino decay also allowed a 
late transfer of energy from non-relativistic to relativistic 
particles, altering the expansion rate history |l4|. How- 
ever, the case considered was that of a heavy (m„ > 10 
eV) neutrino, enough to be the dark matter, decaying 
into massless neutrinos. Such scenarios are no longer 
possible, given laboratory data on the neutrino mass 
scale and mass differences. For the relevant mass range, 
if decays occur, the parent and daughter neutrinos are 
equally relativistic. The possibilities of neutrino anni- 
hilation and/or self- interaction have been considered in 
scenarios in which neutrinos are the dark matter (l5lll6| . 
When neutrinos are a fraction of the dark matter, the 
signatures of this scenario are more subtle, and have not 
been treated elsewhere. 

Annihilation. — The neutrino annihilation rate is 



becomes 



T = (<rv)n cq , 
with the cross section 0, 
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where y/s is the center of mass energy and 1 — 1 — 
4m 2 / /s. In the non-relativistic limit the annihilation rate 
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where we have used (/3 2 ) ~ 3T/ rriv 

For sufficiently large g, the annihilation rate will be 
larger than the expansion rate until the temperature 
drops well below the neutrino mass. Once T„ < m„, 
the neutrino abundance will become exponentially sup- 
pressed, asymptoting to the equilibrium abundance at 
the freezeout temperature, Tf, defined as the tempera- 
ture at which the annihilation rate is equal to the expan- 
sion rate. If Tf is less than of order m v /7, the neutrinos 
will be suppressed from their nominal abundance by a 
factor greater than 100: they will play no role in sub- 
sequent cosmological evolution. The constraint on the 
coupling g is thus obtained by solving T(Tf) = H(Tf) 
and requiring Tf < m v /7. For this estimate, it is suf- 



ficient to set H(T) ~ H [n m {T/T ) 3 + ft 7 (T/T ) 
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where To is the standard photon temperature, 2.73K, and 
f2 7 = 2.47 x 10~ 5 h~ 2 is the ratio of energy density in pho- 
tons to the critical density 18]. Then, we find that as 
long as g > 10 -5 , the annihilation is complete by Tf, 
with only a negligible amount of neutrinos remaining. 

Note that for g > 10~ 5 , the boson will be brought 
into thermal equilibrium before BBN. The energy den- 
sity of a scalar boson is equivalent to 4/ 7 that of a 
neutrino species. Current BBN limits [l9l Eol El) are 
N° s < 3.3 — 4, so an additional boson is still allowable. 
In the case that the electron neutrinos have a large lep- 
ton asymmetry, even N° s = 7 is permitted, provided the 
extra degrees of freedom do not consist of active neu- 
trinos |a E3- Neutrino-majoron interactions may also 
weaken the constraints on large lepton asymmetries |22| . 

Neutrino-boson energy density. — We shall hence- 
forth assume that g > 10~ 5 , so the neutrinos completely 
annihilate into massless bosons. However, there is still 
a small impact on the distribution of matter in the uni- 
verse today. The energy density in the i>~(f) system differs 
from that of the three massless neutrinos of the canon- 
ical standard cosmological model and from a model of 
three massive non-interacting neutrinos. In particular, 
the epoch of matter domination is delayed in the inter- 
acting neutrino scenario outlined above. This delay leads 
to a small suppression of the matter power spectrum on 
small scales. To explain this suppression, we first com- 
pute the evolution of the energy density in the v — <fi 
system and compare it with the conventional scenarios. 

As the neutrinos annihilate, the common temperature 
of the v — 4> fluid does not simply scale as a -1 . Rather, 
it falls less sharply. To track the temperature evolution, 
we can use entropy conservation. The entropy density of 
the v — <b fluid is 
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FIG. 1: Evolution of the energy density as a function of the 
scale factor a. Heavy curves at top are total energy density 
including matter, photons, and neutrinos; light curves at bot- 
tom are energy density in the neutrino sector (including </>'s in 
the interacting case). Three different scenarios are depicted, 
differing in neutrino content: three massless neutrinos (solid), 
three degenerate standard model neutrinos with Y m v = 1 
eV (dotted); and three interacting degenerate neutrinos plus 
massless <j> (dashed). We use the same total matter density, 
Q m =0.3, throughout; p cr denotes the critial debsity today. 
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When the neutrinos are highly relativistic, F = 1, while 
it is exponentially suppressed, F ~ 0, at late times when 
the neutrinos become non-relativistic. Entropy conserva- 
tion then implies 
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If (T^/T-y) i n it takes the standard value, (4/11) 1 / 3 at 
early times, at late times we have (T^^/Tj) = (25/11) 1 / 3 . 
This implies an increase in the radiation energy den- 
sity, corresponding to an effective number of neutrinos 
of Nf = 6.6. The evolution of the energy density is 
shown in Fig. ^ 

CMB measurements constrain the number of light rel- 
ativistic degrees of freedom. The current limit is iV® < 
7 [23, hence does not rule out this scenario. Further, one 
must be careful about applying this limit to our model, 
as interactions will reduce the propagation speed of neu- 
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FIG. 2: The ratio of power spectra P/P(m v = 0) where 
P(m v = 0) is the power spectrum for the standard scenario 
with massless neutrinos. The solid curves show this ratio for 
various (degenerate) neutrino masses in the interacting sce- 
nario. Dashed curves show the ratio in the standard scenario, 
for which the current limit is m v < 1 — 2 eV. Note that the 
tritium bound, m v < 6.6 eV, always applies. 



trinos. Some secondary effects on the CMB due to neu- 
trino freestreaming (i.e., a phase shift and amplitude re- 
duction) will thus be less striking than in the standard, 
non-interacting model |2^]. These effects are discussed 
in |25| for a similar model in which a light (but heavier 
than m v ) boson is coupled to the neutrinos. 

Power spectrum. — We have calculated the large 
scale structure power spectrum, assuming the limit where 
the neutrino annihilation is complete. We find that the 
current neutrino mass limits can be completely removed: 
all values of m v are allowed, even those much greater 
than 1 eV. The result are shown in Fig. [21 where for com- 
parison we have also shown the suppression caused by 
free streaming in the standard case. 

In the interacting scenario, the usual suppression due 
to neutrino mass is absent, because neutrinos make no 
contribution to the matter density today. A small sup- 
pression does occur, due to the extra radiation present. 
Even though neutrinos do not freestream, perturbations 
in the neutrino-0 fluid still cannot grow, due to the pres- 
sure in this tightly coupled relativistic fluid. The negligi- 
ble density in neutrinos makes this suppression irrelevant. 
The effects on the power spectrum are thus entirely due 
to the modified expansion history. 

Matter radiation equality is delayed, since the <f> heat- 
ing leads to an enhanced radiation density (see Fig. 1). 
Therefore, the potentials for scales which enter the hori- 
zon during the radiation dominated epoch will decay for 
a slightly longer period, leading to a small suppression 
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of the power spectrum on these scales. Note that if the 
neutrino annihilation is complete well before matter ra- 
diation equality, as would be the case for very heavy 
neutrinos, the full effects of the extra radiation are felt. 
This corresponds to the bottommost of the solid curves 
in Fig[3 For very small neutrino masses, m„ <C leV, the 
increase in the radiation density due to neutrino annihi- 
lation occurs after time of matter radiation equality. At 
this stage, the universe has already entered the matter 
dominated regime, where the potentials are dominated 
by the dark matter, and the radiation is less important. 
The effects of the extra radiation created by neutrino an- 
nihilation are thus quite small. (The power spectrum is 
slightly suppressed with respect to a standard massless 
neutrino scenario, since there is still a small amount of 
extra radiation due the population of </>.) For interme- 
diate cases, e.g., = leV, we find a suppression 
P/P{m u — 0) ~ 0.8, compared to 0.5 in the normal case. 

Conclusions. — We have examined a model in which 
extra couplings allow the neutrinos to annihilate into 
massless (or light) bosons at late times, and thus make a 
negligible contribution to the matter density today. This 
evades the present neutrino mass limits arising from large 
scale structure. Future tritium beta decay experiments 
like KATRIN will play a unique and essential role, 
especially in comparison to cosmology and neutrinoless 
double beta decay, allowing stringent tests of new neu- 
trino interactions. 

The scenario outlined here could be falsified in several 
ways. First, by a robust discovery of neutrino mass with 
large scale structure data, if the power spectrum suppres- 
sion was greater than that allowed for the tritium bound 
mass in the interacting case (see Fig.EJl. This emphasizes 
the importance of improving the tritium bound. Second, 
with future precision BBN and CMB data. Third, these 
couplings could lead to neutrino decay over astronomical 
distances, which has testable consequences [jfij . 
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